Quasiparticles and c-axis coherent hopping in high T c superconductors 
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We study the problem of the low-energy quasiparticle spectrum of the extended t-J model and 
analyze the coherent hopping between weakly coupled planes described by this model. Starting with 
a two-band model describing the CuO planes and the unoccupied bands associated to the metallic 
atoms located in between the planes, we obtain effective hopping matrix elements describing the 
c-axis charge transfer. A computational study of these processes shows an anomalously large charge 
anisotropy for doping concentrations around and below the optimal doping. 



PACS numbers: 74.25.Fy, 74.20.Mn 

The anomalous charge transport of high T c materi- 
als is still among the most important open problems re- 
lated with the physics of these compounds §],§||. It is 
now accepted that in order to build a complete theory of 
high temperature superconductivity, a deep understand- 
ing of the normal state properties of these compounds 
is needed and one of the most intriguing questions con- 
cerns precisely the charge dynamics ||. The cuprate 
high T c superconductors have crystalline structures con- 
sisting of CuO planes which makes all of them highly 
anisotropic materials. The electrical conductivity reflects 
the anisotropy: while in the plane it shows a metallic be- 
havior, along the c-axis in the underdoped regime, it in- 
dicates an incoherent charge transport; moreover, it has 
an anomalous temperature behavior, and its frequency 
dependence is not of the Drude type Q . 

Band structure calculations together with a semiclas- 
sical theory of charge transport, predict a metallic be- 
havior in all directions with an anisotropy much smaller 
than the experimental values. Electron-electron inter- 
actions may strongly renormalize the bare parameters, 
however within the framework of the Fermi-liquid theory 
the observed anisotropy cannot be understood. This led 
Anderson and others to propose the failure of the con- 
ventional Fermi-liquid theory in these compounds § . 

An important amount of experimental data has been 
accumulated during the last years (|]. Among other 
things, the formation of charge stripes in the normal state 
has been confirmed [Q. Whether this is relevant for the 
occurrence of superconductivity is not yet clear, however 
it should be taken into account in a detailed description 
of the normal phase. 

Here we revisit the problem of the low-energy quasi- 
particle spectrum of the extended t-J model and analyze 
the coherent hopping between weakly coupled planes. 
Our starting point is a stack of CuO planes coupled via 
the electron hopping t±(k) to higher energy intermedi- 
ate states. The intermediate states represent the unoc- 
cupied bands associated to the metallic atoms located in 
between the planes. 

In a single-particle theory, the intermediate states can 



be easily eliminated through a canonical transformation 
giving rise to an effective hopping between consecutive 
CuO planes of the form ^(k) = t 2 ± (\t)/AE k , with AE k 
the one-particle energy difference between k-states of the 
CuO and the intermediate layer bands. The principal 
symmetry of t ± (k) is obtained from electronic structure 
calculations and for the case of YBaCuO it is given by 
[cos(fc x a) — cos(k y a)] 2 ||. 

In what follows we use an uncorrelated empty band 
to describe the intermediate states and the extended t- 
J model to describe the CuO planes. The extended 
t-J model including hole hopping to second and third 
nearest- neighbor (NN) sites, and with J/t consistent with 
experimental values, reproduces the photoemission re- 
sults of the parent insulating materials |^|l0). In the 
standard notation the model Hamiltonian describing the 
plane I is: 

Hi = J ^ Su ■ Sij — ^ t nm cj ncr ci ma (1) 

n.m 

where the matrix element t nm is t for NN, t' for sec- 
ond NN and t" for third NN, q mcr destroys an elec- 
tron with spin a at site m and double occupation is 
not allowed. Using the many-body states, we perform 
a canonical transformation that eliminates the interme- 
diate states and generates an effective Hamiltonian that 
includes a charge transfer matrix element between two 
consecutive CuO planes. Let us first consider the case of 
a single hole. The two-plane states of the system are in- 
dicated as \tpg v ) <g> |'0i v i) where \4> l m „) is the many-body 
wave function of plane I with m holes and quantum num- 
bers v and the corresponding energy is E mM . We take 
the plane with no holes in its ground state indicated by 
v = 0. The effective matrix elements connecting the 
states I^q ) (g> \ip 2 „) and l^i,,,) <8> | Vo,o) accounts for a 
hopping - from one plane to the other - of a single hole. 
The one-hole quantum numbers v include the crystal mo- 
mentum k, spin variables and other quantum numbers u 
that uniquely identify the state and in what follows we 
take v = (k, er, /i). Due to the momentum and spin con- 
servation of the Hamiltonian, the effective hopping of a 



1 



hole with momentum k and spin a is given by: 
t 2 ± (k)M 1 (k, a,fj,)M*(k, (7,/i) 



f ± (k) 



Eq,0 — ^l^kjCT./i) — £ kcr 



(2) 



where £ker is the energy to add an electron at the in- 
termediate band, Mi(k, ct, /i) = (ipl^a^l^^Ofi) and 

M 2 *(k,(r,/x) = (V ? o,ol4k ( rlV'i,(k, -, A1 )> are the matrix ele- 
ment for the creation and destruction of a hole in the 
planes. We calculate many-body energies and matrix el- 
ements by exact diagonalization of small clusters with 
periodic boundary conditions. 

Consider the one-hole ground state in a 16-site clus- 
ter: this state has momentum k = (tt/2,tt/2) and the 
product of the matrix elements in Eq.(|J) is the weight 
of the quasiparticle peak Zk- Since the same factor Zk 
renormalizes the in-plane hole mass, for one hole close 
to the Fermi surface, interactions do not renormalize the 
anisotropy. The one- hole state with k = (7r/2,7r/2) can 
be viewed as a spin-1/2 and charge e quasiparticle, where 
the spins across the hole have ferromagnetic correlations. 
To distinguish this excitation from others that appear at 
different regions of the Brillouin zone (BZ), we refer to it 
as the a-type spin-1/2 quasiparticle. This quasiparticle 
exists for all the momenta in the BZ and its energy can 
be obtained simply by following the lowest energy peak of 
the one-particle spectral density, i.e. its dispersion rela- 
tion, shown in Fig. ([[]), gives the one- hole photoemission 
(PES) dispersion of the insulator. It is important to note 
that this quasiparticle is not the lowest energy state for 
all momenta k. For realistic parameters that fit the PES 
(J/t = 0.4 , t'/t = -0.35 and t" /t = 0.25), the lowest 
energy states in the 4x4 cluster with k = (0, 0) and 
(7r, 7t) have total spin 3/2 p3. These high spin states 
have been obtained in clusters of different sizes (up to 
20 sites) and for different parameters. However it is not 
clear that in the thermodynamic limit they will remain 
being the lowest energy excitations around the T and 
M points. These states can be interpreted as high spin 
quasiparticles that occur as a precursor of the Nagaoka 
state. The magnetization is confined around the hole like 
in a polarized spin bag. These states are not connected 
to the ground state of the undoped system when a hole 
is created; in other words the matrix elements in Eq.(||) 
are zero and these excitations do not propagate - to low- 
est order - from one plane to the other. The reason is 
that when a hole is created on the spin-zero ground state 
of the insulating system, the resulting states have total 
spin 1/2, and the spin wave function is orthogonal to the 
spin-3/2 excitations. 

More interesting are the one-hole lowest energy states 
at k = (tt, 0) or (0, tt). These states have been studied re- 
cently in some detail (n],[ll[] ■ The spin correlations across 
the hole are antiferromagnetic (AFC AH), like in ID sys- 
tems, and it has been speculated that this is an indication 
of charge and spin separation in the 2D t-J model p2[. 



Similar to the high spin excitation case, here also, the 
matrix elements connecting this state with the ground 
state of the undoped system by creating a hole is zero. 
Although the total spin of the k = (tt, 0) lowest energy 
excitation is 1/2, the symmetry of its spin wave function 
is orthogonal to that of cikalV'o.o)- The one-hole spectral 
density of the insulator does not show any structure at 
the energy of this state as shown in Fig. (||d). Again, 
to lowest order, this excitation does not propagate from 
one plane to the other. This is a property of the tt't" — J 
model and it is not observed in the t-J model for similar 
values of J/t. 
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FIG. 1. Dispersion of one hole for the extended t-J model 
on a 4x4 cluster; t'/t = -0.35, t" /t = 0.25, and J/t = 0.4. In 
the inset the one hole PES is shown together with the S — 3/2 
and the AFCAH state energies. Here T = (0,0), M = (7r,7r), 
and X = (tt,0). 

In order to qualitatively determine the regions of the 
BZ where the different type of excitations are the lowest 
energy ones, we have continuously changed the boundary 
conditions in the cluster and the results are shown in Fig. 
(1). In a region around the X point, the spin-1/2 states 
with AFCAH are stable (dashed line). If these excita- 
tions actually describe a region of charge and spin sepa- 
ration, the dashed line branch does not correspond to the 
dispersion of a single quasiparticle but to the superposi- 
tion of a holon and a spinon. It has been suggested that 
these states play a central role in the many-holes states 
and may be responsible for the formation of stripes [p"2| . 
As we show below, they are also relevant for the coherent 
hopping along the c-axis in the doped systems. 

Summarizing the one-hole results, we have shown that 
in different regions of the BZ, the lowest energy exci- 
tations are of different nature. Close to the Fermi sur- 
face (k = (7t/2, 7r/2)), there is a spin-1/2 quasiparticle 
that propagates along the c-axis with an effective hop- 
ping renormalized with the Zk F factor. At the X point, 
the anomalous spin-1/2 excitation with AFCAH is the 
lowest energy state. A hole in this state is confined, up 
to second order in t±, in a single plane. Here the word 
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confinement is used in the sense that there is no bonding- 
antibonding splitting of the degenerate two-plane lowest 
energy state. In Fig.(|]c) the matrix element for different 
values of k is shown. Higher order terms may generate 
a coherent hopping for this type of excitations. In any 
case, we expect these higher order terms to generate a 
very small dispersion along the c-axis jn|. Moreover, 
it should be noted that for the 20-site cluster, for the 
k = (0, 7r) state with robust AFCAH, the second order 
matrix elements non zero but very small, i.e. the quasi- 
particle peak in the spectral density of the insulator has 
an extremely small weight |l0],[ll]], leading to a strongly 
renormalized c-axis hopping for this state, in qualitative 
agreement with the results for the 16-site system. The 
non-crossing approximation in larger systems also gives 
zero weight |l4|] . 

Let us now turn to the many-holes states. For the 
same parameters as above, the two-hole ground state of 
the 4x4 cluster has total momentum k = (0, tt) or (tt, 0) 
and spin zero. In these states, the hole-hole correlation 
function (niJij) with rij the hole number at site i, clearly 
shows the tendency to stripe formation. For the state 
with k = (0, 7r) [(tt, 0)], the holes are mostly aligned in 
the same row along the x [y] direction as shown in Fig 
(2d). These results are in agreement with previous re- 
sults obtained in ladders and clusters. In the 18-site 
cluster, stripes are observed in the three holes ground 
state p^ |. Based on the spin correlations across the 
hole, it has been argued that these stripes are made out 



of one- hole building- blocks with AFCAH |12|. In other 
words, the two holes wave function has a large compo- 
nent of the one-hole wave function of the k = (0, tt) 
type. This could strongly influence the charge mobil- 
ity along the c-axis. We first analyze the case of two 
4x4 planes with three holes corresponding to a doping 
x ~ 0.094. With a straightforward extension of the pre- 
vious treatment and notation, we calculate the effective 
matrix element mixing the two (degenerate) states with 
wavefunctions \^\^ v ) <8> \1p24) anc ^ 1^1 ® l^i,r/) which 
now are Mi(k,z/,£) = (^ e |ci kCT |V>i,„) and M 2 *(k,z/,£) = 



(V'LIcLJV'; 



2kerl V2,{/ 



For these matrix elements to be non 



zero, the operators ciko- and c\^ a have to be chosen to 
conserve total momentum and spin. 

If the two planes are in their ground state, the product 
of these matrix elements is smaller than for the one-hole 
case. This result, is consistent with the conjecture that 
the two holes wave function is mainly made out of one- 
hole building blocks with AFCAH since, as we showed 
above, this type of excitation tends to be confined in one 
plane. This effect is largely enhanced when doping in- 
creases as in the case of five holes in two 4x4 planes 
(x ~ 0.156). In this case, if the planes are in the ground 
state, with two and three holes respectively, the effective 
matrix element to transfer a charge from one plane to 
the other is exactly zero and to lowest order, charge is 



confined in the planes. In fact as shown in Fig. (3b), 
the PES of the system with two holes has no structure 
at the energy of the three-hole ground state. The case of 
four holes (x — 0.125) deserves a special comment: the 
ground state of the uncoupled planes has two holes in 
each one and the hopping mixes it with excited states. 
It has been proposed to artificially shift the energy lev- 
els of the states with one hole in one plane and three 
in the other to recover a degeneracy p5| . These shifts 
are unimportant when calculating the matrix elements 
which are (^l io \c lk(7 \^ o ) x (^i,, l4i«rlV'2,^) wit h the 
indices £o, Vo an d vo corresponding to the ground state 
of a plane with the indicated number of holes. These ma- 
trix elements are precisely the ones calculated in the two 
previous cases and the product is zero. For the case of 6 
holes in two 4x4 planes we reach similar conclusions. 




FIG. 2. (a),(b) Single hole spectral functions A(k,Lj), for 
the extended t-J model, same parameters as Fig.(l) (broaden- 
ing 8 = O.lt). The arrow points to the one-hole lowest energy 
state with the corresponding momentum, (c) Matrix element 
\M\ 2 for the available momenta on a 4 x 4 cluster, (d) Hole 
density around one hole (open circle) on the two holes ground 
state (k = (0,7r)). Full circles radius are proportional to the 
hole density. 

Finally we have evaluated the effective hopping for 
large number of holes, 9 and 11 holes in the two planes 
which correspond to dopings x larger than the optimal 
doping. In this region, the effective hopping calculated 
always in second order in t± and in the ground state, are 
non-zero. 

We have also studied the effective hopping along the c- 
axis for 18-site and 20-site clusters, and small doping. For 
the 18-site cluster, and the same parameters as before, 
the two holes low energy spectrum presents a quaside- 
generacy between the spin zero and spin one states. The 
effective hopping calculated in the lowest energy S = 
state, which we believe is the relevant in the thermody- 
namic limit, shows the same behavior as in the 4x4 
cluster, indicating that our results are not dominated by 
the particular symmetry of these clusters. 
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In Fig. ||(c) we summarize the results for the matrix 
elements, which give an estimate of ratio t±/t° L as a func- 
tion of the doping x, here with is the bare hopping 
obtained from band structure calculations. Our results 
show that while a single hole propagates along the c- 
axis, for doping close and smaller than the optimal dop- 
ing x op ~ 0.2 the lowest order hopping evaluated in the 
ground state cancels |l6| and increases for x > x op . This 
behaviour suggest that when stripes are stabilized, the 
c-axis hopping reduces to very small values. These re- 
sults are in agreement with previous calculations 
which suggest that the weight of the quasiparticle peak 
decreases in the presence of stripes. As we mentioned 
above, higher order terms generate a non-zero hopping. 
However, at least in the one-hole case in the two -\/8x \/8 
planes |fl3|| , these contributions are very small reducing 
the bare hopping by a factor smaller than 10~ 3 . 
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FIG. 3. (a) Single-hole spectral function A(k, u) for a 4 x 4 
cluster with one hole (k = (— 7r/2,7r/2)), for the extended t-J 
model and same parameters as Fig (1). The arrow points to 
the final state energy, (b) Same as (a) for a system with two 
holes (k = (0,0)). (c) Effective hopping between CuO planes 
up to second order in t± as a function of the hole doping x 
for the 16-site and 18-site clusters, (d) Drude weight averaged 
over both directions for the 4x4 cluster as a function of the 
hole doping x. 

For the sake of completeness we have also calculated 
the in-plane Drude weight D for the same parameters as 
before. As shown in Fig. ||(d) D increases with doping 
(as in the t-J model in the region were t± presents a 
marked depletion. Our treatment, based in a canonical 
transformation and the estimation of matrix elements by 
exact diagonalization of small clusters is not appropri- 
ate to give a definite answer to what may be the most 
relevant question related to the charge dynamics in high 
T c materials, that is whether under some conditions real 
confinement of charge excitations in the CuO planes is 
obtained. However our results are consistent with, and 
can explain, the available experimental data. The be- 
havior here obtained is robust under changes of the pa- 



rameters within the accepted range of values. On one 
hand, to lowest order in t±, in the ground state charge 
is confined for doping x close to and below the optimal 
doping x op . Higher order terms in t± may generate a 
small coherent hopping t± ~ t± x 1CP 3 . Estimations of 
the bare hopping using LDA band structure calculations 
give t± = 0.25 eV for the case of bilayers and t± ~ 0.03 
cV for the intcrplane matrix element B. That means 
that down to very low temperatures (kT ~ t±), we ex- 
pect a diffusive-like propagation along the c-axis with an 
anomalous temperature dependence. On the other hand, 
the small cluster calculations are remarkable in the sense 
that they predict an anomalously large charge anisotropy 
precisely in the region where superconductivity is stable. 
This agrees with the experimental data and strength- 
ens the ideas developed by Anderson and coworkers [HJ 
that argued that the lack of single particle intraplane co- 
herent hopping enhances superconductivity. 
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